home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlags2.z / dlags2
Encoding:
Text File  |  2002-10-03  |  4.1 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAGGGGSSSS2222((((3333SSSS))))                                                          DDDDLLLLAAAAGGGGSSSS2222((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAGS2 - compute 2-by-2 orthogonal matrices U, V and Q, such that if (
  10.      UPPER ) then  U'*A*Q = U'*( A1 A2 )*Q = ( x 0 ) ( 0 A3 ) ( x x ) and
  11.      V'*B*Q = V'*( B1 B2 )*Q = ( x 0 ) ( 0 B3 ) ( x x )  or if ( .NOT.UPPER )
  12.      then  U'*A*Q = U'*( A1 0 )*Q = ( x x ) ( A2 A3 ) ( 0 x ) and V'*B*Q =
  13.      V'*( B1 0 )*Q = ( x x ) ( B2 B3 ) ( 0 x )  The rows of the transformed A
  14.      and B are parallel, where  U = ( CSU SNU ), V = ( CSV SNV ), Q = ( CSQ
  15.      SNQ ) ( -SNU CSU ) ( -SNV CSV ) ( -SNQ CSQ )  Z' denotes the transpose of
  16.      Z
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  19.      SUBROUTINE DLAGS2( UPPER, A1, A2, A3, B1, B2, B3, CSU, SNU, CSV, SNV,
  20.                         CSQ, SNQ )
  21.  
  22.          LOGICAL        UPPER
  23.  
  24.          DOUBLE         PRECISION A1, A2, A3, B1, B2, B3, CSQ, CSU, CSV, SNQ,
  25.                         SNU, SNV
  26.  
  27. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  28.      These routines are part of the SCSL Scientific Library and can be loaded
  29.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  30.      directs the linker to use the multi-processor version of the library.
  31.  
  32.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  33.      4 bytes (32 bits). Another version of SCSL is available in which integers
  34.      are 8 bytes (64 bits).  This version allows the user access to larger
  35.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  36.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  37.      only one of the two versions; 4-byte integer and 8-byte integer library
  38.      calls cannot be mixed.
  39.  
  40. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  41.      DLAGS2 computes 2-by-2 orthogonal matrices U, V and Q, such that if (
  42.      UPPER ) then U'*A*Q = U'*( A1 A2 )*Q = ( x 0 ) ( 0 A3 ) ( x x ) and
  43.      V'*B*Q = V'*( B1 B2 )*Q = ( x 0 ) ( 0 B3 ) ( x x ) or if ( .NOT.UPPER )
  44.      then U'*A*Q = U'*( A1 0 )*Q = ( x x ) ( A2 A3 ) ( 0 x ) and V'*B*Q = V'*(
  45.      B1 0 )*Q = ( x x ) ( B2 B3 ) ( 0 x ) The rows of the transformed A and B
  46.      are parallel, where U = ( CSU SNU ), V = ( CSV SNV ), Q = ( CSQ SNQ ) (
  47.      -SNU CSU ) ( -SNV CSV ) ( -SNQ CSQ ) Z' denotes the transpose of Z.
  48.  
  49.  
  50. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  51.      UPPER   (input) LOGICAL
  52.              = .TRUE.: the input matrices A and B are upper triangular.
  53.              = .FALSE.: the input matrices A and B are lower triangular.
  54.  
  55.      A1      (input) DOUBLE PRECISION
  56.              A2      (input) DOUBLE PRECISION A3      (input) DOUBLE PRECISION
  57.              On entry, A1, A2 and A3 are elements of the input 2-by-2 upper
  58.              (lower) triangular matrix A.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAAGGGGSSSS2222((((3333SSSS))))                                                          DDDDLLLLAAAAGGGGSSSS2222((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      B1      (input) DOUBLE PRECISION
  75.              B2      (input) DOUBLE PRECISION B3      (input) DOUBLE PRECISION
  76.              On entry, B1, B2 and B3 are elements of the input 2-by-2 upper
  77.              (lower) triangular matrix B.
  78.  
  79.      CSU     (output) DOUBLE PRECISION
  80.              SNU     (output) DOUBLE PRECISION The desired orthogonal matrix
  81.              U.
  82.  
  83.      CSV     (output) DOUBLE PRECISION
  84.              SNV     (output) DOUBLE PRECISION The desired orthogonal matrix
  85.              V.
  86.  
  87.      CSQ     (output) DOUBLE PRECISION
  88.              SNQ     (output) DOUBLE PRECISION The desired orthogonal matrix
  89.              Q.
  90.  
  91. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  92.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  93.  
  94.      This man page is available only online.
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.